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Abstract 

We discuss inclusive production of open charm in proton-proton scattering at LHC. The calcu- 
lation is performed within the ^-factorization approach. Different models of unintegrated gluon 
distributions (UGDF) from the literature are used. The theoretical transverse momentum as well 
as (pseudo) rapidity distributions of charmed mesons are compared with recent experimental data 
of ATLAS, ALICE and LHCb collaborations. Only the calculation with Kimber-Martin-Ryskin 
(KMR) UGDF gives results comparable to experimental ones. All other popular models of UGDF 
significantly underpredict experimental data. Several sources of uncertainties of the theoretical pre- 
dictions are also studied in details. In addition we discuss correlations between D and D mesons. 
Good description of experimental distribution in invariant mass and in relative azimuthal angle 
between D and D mesons is achieved for the KMR UGDF. The considered correlation observables 
measured by the LHCb experiment were not discussed in other approaches in the literature. 
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I. INTRODUCTION 



At high energy hadronic scattering the gluon-gluon fusion is known to be the dominant 
mechanism of open charm production. Even at RHIC the contribution from quark- ant iquark 
anihilation constitutes only a small fraction of the cross section. Usually in the studies of 
heavy quark production the main efforts concentrate on inclusive distributions. The trans- 
verse momentum distribution of charmed mesons is the best example. Standard collinear 
NLO approach [l[ as well as its improved schemes e.g. FONLL Q or GM-VFNS 0] are states 
of art in this respect. These approaches cannot be, however, used when transverse momenta 
of charm quark and antiquark are not equal. This means in practice that it cannot be used 
for studies of correlation observables for charmed meson pairs or for meson-nonphotonic 
electron modes. 

The k t - factorization approach seems much more efficient tool in this respect HJ- Differ- 
ent unintegratedgluon distributions in the proton (UGDF) have been used in the literature 
in this context 12-16||. Recently we have applied this formalism to the description of inclu- 



sive distributions of so-called nonphotonic electrons [17j and electron-positron correlations 



181 ] at RHIC. Rather good description of correlation observables has been achieved there. 



The quark mass is sufficiently large to apply perturbative calculation, but still small 
enough that interesting low-x effects may appear too. In the fc t -factorization approach the 
latter effects are contained in the unintegrated gluon distributions - the buildning blocks 
of the formalism. In principle a comparison of experimental data and predictions with the 
UGDFs which include such effects may tell us more about footprints of the saturation effects 
- the topic extensively discussed in recent years. 

Recently ATLAS 0, ALICE HQ and LHCb Q collaborations have measured inclu- 
sive distributions (mainly transverse momentum distributions) of different charmed mesons. 
The LHCb collaboration has measured in addition a few correlation observables for charmed 
mesons for the first time in the history in the forward rapidity region [23|]. Before the STAR 
collaboration at RHIC has measured correlation of charmed mesons and nonphotonic elec- 



trons |24|. At RHIC a study of meson-meson correlations was not possible due to limited 
statistics caused by relatively small cross sections. It was accessible only at Tevatron where 
first midrapidity measurements of azimuthal ang le correlations between charmed mesons 



have been performed by the CDF experiment [25 . 

In the present paper we wish to concentrate first on inclusive distributions of charmed 
mesons in order to test different models of unintegrated gluon distributions from the litera- 
ture. Next we wish to focus on DD meson correlations. Conclusions will close our paper. 



II. SKETCH OF THE FORMALISM 

The cross section for the production of a pair of charm quark - charm antiquark can be 
written as: 

dajpp -> ccX) 1 f d 2 k lt d 2 fa t , .. of f ^ 

d yi dy 2 d 2 p lt d 2 p 2t ~ = 16^P J ~ ~ l 9 * 9 ^ cA 

x 5 2 (ku + fat -Pit ~P2t) J r g (xi,kl t ,fi 2 )J : ' g (x2,kl t ,iJ 2 ). (2.1) 

The main ingredients in the formula are off-shell matrix elements for g*g* — y c c subprocess 
and unintegrated gluon distributions (UGDF). The relevent matrix elements are known 
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and can be found in Refs. 126-281. The unintegrated gluon distributions are functions of 
longitudinal momentum fraction x\ or x 2 of gluon with respect to its parent nucleon and of 
gluon transverse momenta k t . Some of them depend in addition on the factorization scale 
fi. The longitudinal momentum fractions can be calculated as: 

xx = -=exp(j/i) + -^exp(y 2 ), 
V s V s 

X2 = —/= exp(-yi) + — = exp(-y 2 ), (2.2) 



where m it = yPu + m Q is the transverse mass of produced quark/ ant iquark. 

Various unintegrated gluon distributions have been discussed in the literature [12|-|l6| . In 
contrast to the collinear gluon distributions (PDFs) they differ considerably among them- 
selves. One may expect that they will lead to different production rates of cc pairs at the 
LHC Since the production of charm quarks is known to be dominated by the gluon-gluon 
fusion, the charm production at the LHC can be used to verify the quite different models of 
UGDFs. 

Below we wish to concentrate for a while on the Kimber-Martin-Ryskin (KMR) uninte- 
grated gluon distribution which, as will be discussed in this paper, gives the best description 
of the LHC experimental data, taking into account also correlation observables. 

According to the KMR approach the unintegrated gluon distribution is given by the 
following formula 

f g (x } k 2 } ^ 2 ) = - ^ [ g ( Xj k 2 t ) T g (k 2 , /1 2 ) ] 

= T g (kl^)^J2jj zP ^ b {- z ^- ( 2 - 3 ) 

This definition is fully satisfied for k t > /io, where fio ~ 1 GeV is the minimum scale for 
which DGLAP evolution of the conventional collinear gluon distributions, g(x,fi 2 ), is valid. 

The virtual (loop) contributions may be resummed to all orders by the Sudakov form 
factor, 



T g (k 2 ^ 2 ) = exp ( - £ M ^Ml ^ / dzzPbg{z) ) , (2.4) 



J2 J o dz z P bg (z)j , 



which gives the probability of evolving from a scale k t to a scale fi without parton emission. 

The exponent of the gluon Sudakov form factor can be simplified using the following 
identity: P qg (l — z) = P qg (z). Then the gluon Sudakov form factor is 



T g (k 2 ^ 2 ) = exp (-£ ^ ^1 ^ A dzzP gg (z) +n F j\zP qg {z)^ 



(2.5) 



where np is the quark-antiquark active number of flavours into which the gluon may split 
and A = k t /(k t -\- ti) which introduces a restricton of the phase space for gluon emmision due 
to the angular-ordering condition. Due to the presence of the Sudakov form factor in the 
KMR prescription only last emission generates transverse momentum of incoming gluons. 
This scheme is the direct analogy to the techniques usually applied in all standard parton 
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shower Monte Carlo generators. The unique feature of the KMR model of UGDF is that 
it provides possibility for the emission of at most one additional gluon. Therefore one can 
expect that the KMR model may include in an effective way NLO corrections to heavy quark 
production cross section. 

In the literature often somewhat differently defined UGDFs are used. They differ by the 
following transformation: 



1 



T g (x, k t , fi 2 ) — ^ fg{ x i k 2 , fj 2 ) 
The normalisation condition for unintegrated distributions 



(2.6) 



d^t fg( x i k t , fJ> ) 



(2.7) 



is exactly satisfied if we define 

1 



2 fg{ x i k 2 , fJ 2 ) 



k t </i 



/'5 



g(x,fi 2 Q )T g ({i 2 ,{i 2 



so that the density of gluons in the proton is constant for k t < fi at fixed x and fi. 
The precise expression for the unintegrated gluon distribution reads 



fg(x, k 2 , fj 2 



2tt 



X 



dz 



J2p 9q ( z )- zq (^k 2 ) + p gg ( z ^ 9 (^)e 



ii 



(2.8) 



(2.9) 



III. CHARM QUARK/ANTI-QUARK PRODUCTION AT LHC 

In this section we wish to concentrate on the production of charm quarks and antiquarks. 
Thus this section has rather theoretical character. The cross sections for production of 
charmed mesons will be discussed in the next section. Before we go to the presentation of 
differential distributions let us summarize integrated cross sections for cc production. 

Using the KMR model of unintegrated gluon distributions, the total cross section for 
charm quark/antiquark production at yfs = 7 TeV is obtained to be cr^ IR {pp — > ccX) = 
7.36^)1* '%l(m c ) mb. The predicted value has large uncertainties related to the choice of 
factorization/renormalization scales /i and due to the charm quark mass m c . The obtained 
cross section is very large, of the same order as e.g. cross section for elastic scattering or 
single diffraction. This means that in practice charm quark/antiquarks appear in almost 
each inelastic event. This is a rather new situation which requires more detailed studies. 

Taking into account acceptance of ATLAS, LHCb and ALICE detectors we get 
<¥l R as(pP ccX) = 2.53l°;^)t^K) mb, -> ccX) = lM^M^m.) 

mb and cr^lj^ E (pp — > ccX) = O-d^o^i^-o^li'rnc) m b, respectively. These numbers to- 
gether with theoretical uncertainties are consistent with recent LHC measurements as well 
as with the recent FONLL [2^] and GM-VFNS j30| predictions of charm cross section. 

As it was mentioned in the previous section our predictions are very sensitive to the choice 
of unintegrated gluon distributions. Different UGDFs are very often based on quite different 
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theoreticl assumptions. This has a crucial meaning for their kinematical characteristics. 
In Fig. [I] we show dependence of the unintegrated gluon distributions functions on gluon 
transverse momentum squared kf for several values of x relevant for the production of charm 
quarks and antiquarks at LHC energy. Differences in shapes in kf of the plotted functions are 
significant. One can also see different dependence on x of the different considered UGDFs. 
Changing the value of x, the mutual trends between them also change what makes the 
overall picture more complicated. Especially the KMR model seems to reveal the strongest 
x- dependence. 




FIG. 1: Different unintegrated gluon distributions from the literature as a function of gluon 
transverse momentum squared kj for different values of longitudinal momentum fraction x of the 
gluon initiating the hard process and for different factorization scale [i. 



The rapidity of the quark or antiquark are strongly correlated with longitudinal momen- 
tum fractions of gluons initiating the hard process. This is shown in Fig. [2] for the KMR 
UGDF. At rapidities \y\ > 5 one starts to probe longitudinal momentum fractions smaller 
than 10 -4 . This is a new situation compared to earlier measurements at RHIC or Tevatron. 
The unintegrated gluon distributions (UGDFs) as well as standard collinear ones (PDFs) 
were not tested so far in this region. 



p p -> c c X is = 7 TeV p p -> c c X is = 7 TeV 




log (x ) log (x ) 

FIG. 2: The range of longitudinal momentum fraction of gluons and its correlation to the rapidity 
of charm quark (left) or antiquark (right). In addition regions of the coverage for the ATLAS and 
LHCb experiments are shown. 
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It was advocated in Ref. 31( that the two-dimensional distribution in transverse mo- 
mentum of charm quark and charm antiquark can be a good "theoretical observable" to 
study unintegrated gluon distributions. In Fig. [3] we show such distributions for different 
UGDFs from the literature. We use here KMR Q, KMS 0, Kutak-Stasto Q, Jung 
setA+, setB+ [15] and GBW [l6| parametrizations. Quite different pattern is obtained for 
different UGDFs. This may have direct consequences for correlation observables for mesons 
or/and nonphotonic electrons. Moreover, events when one p t is small and second one is 
large correspond to the region relevant for higher order collinear corrections. It is clear from 
this pitp2rplane that effects of an effective inclusion of NLO diagrams in the fct-factorization 
approach strongly depend on the construction of UGDFs. 



p p -» c c X 
> (lyjsB '* 

<» OR 

(5 25 ( 



p p -> c c X 



is = 7 TeV 



p p -> c c X 



is = 7 TeV 
Jung setA-H^ 




charm quark (GeV) 



charm quark (GeV) 



charm quark p^ (GeV) 



p p -» c c X 



is = 7 TeV 
■'jffing sctB+ 



p p -» c c X 




p p -> c c X 



is = 7 TeV 



> 

CD „ 

a 2 ~ 



to 



charm quark p^ (GeV) 



charm quark p^ (GeV) 



^^^^^^^ 



charm quark p^ (GeV) 



3 



4 



10 3 
10 2 
10 

1 

10" 1 
10 2 
10" 3 
10 4 



FIG. 3: Two dimensional maps in transverse momentum of charm quark and transverse momen- 
tum of the charm antiquark for different unintegrated gluon distributions. 



The production of charmed mesons strongly depends on the choice of UGDF model as 
will be discussed in the next section. As will become clear there the KMR UGDF within 
rather large theoretical uncertainties provides the best description of the LHC experimen- 
tal data. The major part of these uncertainties comes from the perturbative part of the 
calculation. Therefore in the following we wish to spend some time to define uncertainties 
of the corresponding calculations at the quark level. In Fig. H] and Fig. [5] we present the 
uncertainties of our predictions, obtained by changing charm quark mass m c = 1.5 ± 0.3 
GeV and by varying renormalization and factorization scales fi 2 = (m 2 , where ( G (0.5; 2). 
The gray shaded bands represent these both sources of uncertainties summed in quadrature. 
The smaller transverse momentum the larger uncertainty. For comparison we show also 
results for the FONLL 0] and MC@NLO (denoted on figures as NLO PM) Q approaches. 
Our result of the /c r factorization approach is consistent within the uncertainty bands with 
those rather standard NLO collinear calculations. Only at small quark p t 's some difference 
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appears. This is the region where transverse momenta of incident gluons play an important 
role. Particularly, a detailed treatment of the nonperturbative kt region in UGDF may lead 
to a dumping or an enhancement of the cross section at small p t . 




FIG. 4: Theoretical uncertainties on transverse momentum distribution of c or c production due 
to the choice of factorization/renormalization scale and those related to charm quark mass for 
the KMR UGDF (solid line with the shaded bands). Left panel shows the cross section for the 
whole range of quark/antiquark rapidities while the left panel for the rapidity range relevant for 
the LHCb experiment. For comparison the FONLL and NLO PM predictions are also shown. 




5 10 15 20 25 30 5 10 15 20 25 30 



P ± (GeV) p ± (GeV) 

FIG. 5: The same as in Figfjbut for the ALICE (left panel) and ATLAS or CMS (right panel) 
kinematics. 



IV. PRODUCTION OF CHARMED MESONS 



The hadronization of heavy quarks is usually done with the help of fragmentation func- 
tions. The inclusive distributions of charmed mesons can be obtained through a convolution 
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of inclusive distributions of charm quarks/antiquarks and c — > D fragmentation functions: 



da(pp -> DDX) f 1 dz da(pp ^ ccX) 



f dz < 

Jo z 



dy D d 2 p t ,D Jo z 2 dy c d 2 p 



Jt,c 



(4.1) 



Pt,c=Pt,D/ z 



where p tyC = and z is the fraction of longitudinal momentum of heavy quark carried 
by meson. We have made typical approximation assuming that y c is unchanged in the 
fragmentation process, i.e. Ud = Vc- 

As a default set in our calculations we use standard Peterson model of fragmentation 
function [33) with the parameter e c = 0.05. This value was extracted by ZEUS and HI anal- 
yses and seems to be relevant for LO calculations. Hoowever, in the fragmentation scheme 



applied in the FONLL framework, rather harder functions (or smaller e c ) are suggested [34 



This issue together with effects of applying other fragmentation functions from the literature 



[35H37J will be discussed in more detail when discussing differential distributions. 

In Table [I] we have collected integrated cross sections for the production of different 
species of D mesons. Measured cross sections from different LHC experiments are compared 
to theoretical predictions obtained with three sets of UGDFs. The error bars shown for the 
KMR UGDF reflect uncertainty due to the choice of factorization/renormalization scale (/x) 
and related to the mass of the quark (m c ). The fractional uncertainties due to these both 
sources for other UGDFs are similar. Only cross sections obtained with the KMR UGDF 
are consistent within error bars with the experimental data. 

In the cases of measurements with the full coverage of the meson transverse momentum 
range, the theoretical cross sections are almost insensitive to the fragmentation model. Quite 
different situation is observed when small pt region is excluded. In the latter case, using the 
Peterson model with e c = 0.02 (which gives results closer to the FONLL predictions) we 
note the enhancement of the integrated cross sections by about 20%. 

Let us start presentation of differential distributions for different LHC experiments. 



ALICE 



Let us focus first on the production of charmed mesons at midrapidities. The ALICE 
collaboration has performed a measurement of transverse momentum distribution of D°, D + , 
In the very limited range of (pseudo) rapidity one tests unintegrated gluon 
distributions in a pretty narrow region of longitudinal momentum fractions (see Fig. [2]). In 
Fig. [6] we show transverse momentum distribution of D° mesons. In the left panel we present 
results for different UGDF known from the literature. Most of the existing distributions fail 
to describe the ALICE data. The KMR UGDF provides the best description of the measured 
distributions. Therefore in the following we shall concentrate on the results obtained with the 
KMR UGDF. In the right panel we show uncertainties due to the choice of usual integrated 
collinear gluon distributions (PDFs) used for calculating the KMR UGDF. In the latter 
case the biggest uncertainty can be observed at small transverse momenta, i.e. in the region 
of small gluon longitudinal momentum fraction. We use rather u p-to -date MSTW08 [38], 
CTEQ6 [39] and GJR08 40 1 parametrizations as well as GRV94 41] which is fairly older 
but was very often used in last years in similar analyzes. For more detailed discussion of 
the PDFs aspects in charm production we refer the reader to Ref . |42) . In Fig. [7] we show 
separately uncertainties due to the choice of factorization/renormalization scale (left panel) 



TABLE I: Integrated cross sections for production of different D mesons at LHC. 



Acceptance 


Mode 




KMR t(m c ) 


Jung setA0+ 


KMS 


ALICE 


(D° + D°)/2 


516 ± 41±?? B 


K-IA + 169 + 384 

,J±4: -130 -198 


317 


313 


\y\ < 0.5 


{D+ + D-)/2 


248 ± 30l| 


9DR +68 +154 
zuo -52 -79 


127 


125 




(D* + + D*~)/2 


247 ± 27±|f 


zuo -53 -80 


129 


127 


ALICE 












\y\ < 0.5 


(D+ + D s )/2 


53 ± 12lig 




13 (.+ /u/o i 


13 \ +AU/o ) 


2 < p ± < 12 GeV 












LHCb 


D° + D° 


1488 ± 182 


1 7AA + 565 + 1435 
i'** -418 -700 


1162 


872 


2 < y < 4.5 


T~\-\- i T~\ — 

+ D 


717 ± 109 


f>Q7 +226 +574 
DJ ' -167 -280 


465 


349 


< p ± < 8 GeV 


D*+ + D*- 


676 ± 137 


7nc - +229 +582 
' uo -169 -284 


471 


354 




D+ + D 7s 


194 ± 38 


9/1 fi + 80 + 203 
zqx> -59 -99 


164 


123 


ATLAS 

M < 2- 1 
p± > 3.5 GeV 


D + + D- 

D+ + D s 
D *+ + D *- 


238 ± 13±f| 
168 ± 34+J| 
285 ± 16±p 


i o 7 +31 +30 (+20%) 
10 ' -20 -24 
4 o +12 +11 (+20%) 
^° -7 -8 
ice +37 +37 (+20%) 
100 -22 -28 


103 (+ 20% ) 
36 (+ 2 o%) 

115 (+20%) 


93 (+20%) 
33 (+20%) 

1Q4 (+20%) 



and those due to the choice of the quark mass (right panel). The uncertainties due to the 
choice of scales is rather large. The uncertainties due to quark mass are significant only at 
small transverse momenta. They are calculated by varying the quark mass m c = 1.5 ± 0.3 
GeVand are representative for all other UGDFs. 




is = 7 TeV 

• ALICE 

KMR 

Jung setA+ 

Jung setAO 

Jung setA- 

- Jung setB+ 

KMS 

Kutak-Stasto 



P ± (GeV) 



7 TeV 



ALICE 

— CTEQ6 LO 



GRV94LO 

GJR08 LO 
M ST W0 8 LO 




p ± (GeV) 



FIG. 6: Transverse momentum distribution of D mesons for the ALICE measurement. The left 
panel shows results for different UGDFs while the right panel shows uncertainties due to the choice 
of collinear gluon distributions in calculation of the KMR UGDF. 

In Fig. [8] and Fig. [9] we present corresponding plots for D + mesons. The situation here 
is very similar to the case of D° mesons. 

Let us quantify now uncertainties due to the fragmentation process. Fig. [TOl shows results 
for D° (left panel) and D + (right panel) mesons for different fragmentation functions from 
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FIG. 7: Uncertainties of the theoretical cross section for the D° meson production within the 
ALICE acceptance due to the choice of the scale (left) and due to the quark mass (right). 



pp^D + X 




s = 7 TeV 

ALICE 
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— Jung setA+ 
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>-> Jung setA- 

- Jung setB+ 
KMS 

- Kutak-Stasto 




P_L (GeV) 



P ± (GeV) 



FIG. 8: The same as in Fig. [6] but for the production of D + mesons. 




FIG. 9: The same as in Fig. [7] but for the D + meson. 



the literature. We use here the Peterson model with three different sets of e c parameter, as 
well as Braaten et al. 35|. Kartvelishvili et al. 36| and Collins-Spiller 37J parametrizations. 
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All of the applied functions give similar results. The effects related to the fragmentation 
process seem to be important only at larger meson pt's, starting from p t = 3 GeV. 
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FIG. 10: Uncertainties due to fragmentation parameter e c in the Peterson fragmentation function 
and related to the choice of other fragmentation models for the KMR UGDFs. 

Let us compare now results of our approach to the results of some other popular ap- 
proaches used in the literature. In Fig. CD] and Fig- LEI we present such a comparison. Our 
results obtained within the /c r factorization approach with the KMR UGDF are very similar 
to those obtained within NLO PM and FONLL models. The cross sections obtained within 
leading-order collinear approximation (LO PM) are much smaller, in particular for larger 
transverse momenta. Comparing left and right panels of these figures one can observe an 
improvement of the large pt data description when e c = 0.02 in the Peterson function is 
taken. This choice corresponds to the upper limit of our uncertainties in the hadronization. 
It makes our results closer to those from FONLL. In the FONLL approach as a default 
fragmentation scheme for charm quarks the BCFY model with r c = 0.1 is used. However, 
the Peterson parametrization with e c = 0.02 gives in general very similar characteristics. 
Since our fct-factorization calculation is very similar to the FONLL predictions at the quark 
level (see Fig. [5]) application of the harder fragmentation functions may be justified. 

Now let us consider for the moment distributions for vector mesons D* + . In Fig. [T31 we 
show transverse momentum distributions of D* + for different UGDFs and uncertainties in 
calculating distributions with the KMR UGDF due to the choice of collinear gluon distri- 
butions. As in the previous cases, the choice of collinear PDFs has some importance only 
at small transverse momenta. The same conlcusions as in the cases of pseudoscalar mesons 
come from Fig. dU where uncertainties due to the scales (left) and related to the quark 
mass (right) are presented. In Fig. [TS]both of these sources are taken together and our pre- 
dictions with the KMR UGDF are confronted once again with LO and NLO PM collinear 
calculations. Here we use the Braaten et al. model for fragmentation which has, the only 
one on the market, parametrization of fragmentation function for the transition of heavy 
quark into vector meson state. 

Finally in Fig. [16] we show distributions for mesons, i.e. mesons built of charm and 
strange quarks/antiquarks. The corresponding cross section is considerably smaller than for 
the charm mesons containing light (up or down) quarks/antiquarks. The general situation 
is, however, very similar. The KMR UGDF provides the best agreement with the ALICE 
data. Results for other UGDFs are much below the experimental data which means, in our 
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FIG. 11: Transverse momentum distribution of D mesons for the ALICE kinematical region for 
different values of the Peterson e c parameter. Together with our predictions for the KMR UGDF 
(solid line with shaded band) results of different popular approaches are also shown. 




FIG. 12: The same as in Fig. [IT] but for the D + meson. 



opinion, that they do not pass the powerfull test. In the case of Dg meson a different, quite 
important source of uncertainties appears, namely the poorly known fragmentation fractions 
BR(c -> D+). Changing the value from 0.08 (ZEUS, HI) to 0.116 (ALEPH) a significant 
enhancement of the theoretical predictions is achieved. 



B. ATLAS 



The ATLAS experiment covers much broader range of pseudorapidities than ALICE. As 
a consequence one tests broader region of longitudinal momentum fractions 10 -4 < Xi, x% < 
10~ 2 . The gluon distributions in this range of X\ and x 2 values carried by gluons are rather 
well known. Also application of the known UGDFs from the literature should be reliable. 

Fig. [T7] shows transverse momentum distributions of charged pseudoscalar D ± mesons 
for different models of unintegrated distributions, for e c = 0.05 (left) and e c = 0.02 (right). 
General situation is very similar as for the ALICE experiment although the agreement is 
somewhat worse. Only the upper limit of the KMR result is compatible with the ATLAS 
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P ± (GeV) 



P 1 (GeV) 



FIG. 13: Transverse momentum distribution of D* + mesons for different UGDFs (left) and the 
uncertainties due to the choice of collinear PDFs used in calculating the KMR UGDF (right). 




FIG. 14: Transverse momentum distribution of D* + mesons. Shown are uncertainties due to the 
choice of scales in the KMR UGDF (left panel) and due to the quark masses (right panel). 



experimental data. This may be caused by much broader range of pseudorapidities in the 
case of the ATLAS detector. Potentially, this can be related to double-parton scattering 
effects to be discussed elsewhere {43[. Also the other standard approaches give results below 
the ATLAS data as can be seen in Fig. [181 

Due to fairly large span of pseudorapidities the ATLAS collaboration can extract also 
pseudorapidity distributions. We wish to show now also results for charm meson pseudora- 
pidity distributions. In Fig. [19] and Fig. [20] we show pseudorapidity distributions for charged 
D ± meson. These distributions are rather flat. The results are also compared to recent (pre- 
liminary) ATLAS data. Only the upper limits of large error bars of the theoretical results 
obtained with the KMR distributions are consistent with the ATLAS data. The results with 
other UGDFs clearly underpredict the experimental data. 

As for pseudoscalar mesons above, in Fig. [2T]we show transverse momentum distributions 
for charged vector mesons. The situation is pretty much the same as for pseudoscalar charged 
mesons discussed previously in Fis. [17] and Fig. [18] 

In Fig. [20] we compare results obtained within the /^-factorization approach (grey band) 
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P ± (GeV) 

FIG. 15: Transverse momentum distribution of D* + mesons for different approaches known from 
the literature. The grey band represents overall uncertainties of the /^-factorization approach with 
the KMR UGDF. 




P ± (GeV) p ± (GeV) 



FIG. 16: Transverse momentum distribution of Df mesons. We show uncertainties due to the 
choice of scales and masses (left panel) and uncertainties related to poorly known fragmentation 
branching fraction BR(c — > D$) (right panel). 



with results obtained within other approaches. The central value of the fc r factorization 
approach with the KMR UGDF is consistent with the FONLL and NLO PM predictions. 



C. LHCb 

Finally let us focus on the measurements in the forward rapidity region 2 < y < 4.5. 
Recently the LHCb collaboration presented first results for the production of D°, D + , D* + 



and Df mesons [22j. In this region of phase space one tests asymmetric gluon longitudinal 
momentum fractions: X\ ~ 10~ 5 and Xi > 10~ 2 (see Fig. [2]). This is certainly more difficult 
region for reliable calculation and interpretation of experimental data. First of all gluon 
distributions were never tested at such small x\ values. Secondly many UGDFs from the 
literature may be not good enough for X2 > 10~ 2 . Therefore some care in interpreting the 
results is required. 
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FIG. 17: Transverse momentum distribution of D mesons for different UGDFs from the literature 
compared with the preliminary ATLAS experimental data for e c = 0.05 (left) and e c = 0.02 (right). 




FIG. 18: Transverse momentum distribution of D ± mesons for different standard approaches 
compared with the ATLAS experimental data 19|] for e c = 0.05 (left) and e c = 0.02 (right). The 
details are specified in the figures. 



The LHCb, similar as ALICE, has measured also distributions of rather rarely produced 
.D+ mesons. We start from transverse momentum distributions for Df mesons. In Fig. [221 we 
present distributions for different UGDF from the literature and uncertainties for the KMR 
UGDF related to the choice of standard PDFs. In Fig. [23] we show uncertainties related to 
the choice of factorization/renormalization scale and due to the choice of quark masses and 
in Fig. [2H uncertainties related to fragmentation functions. All these uncertainties are very 
similar as for the ALICE and ATLAS kinematics. 

In Fig. [25] we compare our predictions, together with predictions of other popular ap- 
proaches. Main conclusions are the same again as for the ALICE and ATLAS. 

The LHCb collaboration was able to measure transverse momentum distributions of 
mesons in many narrow bins of (pseudo) rapidity. Below (Figs. [26], [27] and 128]) we show 
such distributions for D°, D + and D* + , respectively. In general, different bins are sensitive 
to different regions of longitudinal momentum fractions carried by gluons. However, we do 
not observe any interesting trend in the quality of the description of the LHCb data. Our 
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FIG. 19: Distribution in meson pseudorapidity. The results for different UGDFs are compared 
with the ATLAS preliminary data 19| for different values of the parameter e c of the Peterson 
fragmentation function: e c = 0.05 (left), e c = 0.02 (right). 
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FIG. 20: Distribution in D ± meson pseudorapidity. The results of our calculation are compared 
with those for other calculations and with the ATLAS preliminary data 19j for different values of 
the parameter e c of the Peterson fragmentation function: e c = 0.05 (left), e c = 0.02 (right). 



results with the KMR UGDF within uncertainties are consistent with the experimental data 
and with the FONLL and NLO PM predictions. Corresponding distributions for PYTHIA 
are taken from Ref. j22| and have slightly different p r slope than the other ones. 
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FIG. 21: Distribution in D* + meson transverse momentum. The results of our calculation are 
compared with the ATLAS preliminary data [lQf j . In the left panel we show results for different 
UGDFs and in the right panel our results are compared with other approaches. 
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FIG. 22: Transverse momentum distribution of Df mesons for different UGDFs from the literature 
(left) and the dependence on the choice of collinear gluon distribution functions for the KMR UGDF 
(right). The results of calculation are compared with the LHCb collaboration data. 
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FIG. 23: Uncertainties of the theoretical predictions due to the choice of scales for the KMR 
UGDF (left) and due to charm quark mass (right). 



pp^(d; + d s )x 

.2 <y n 



s 

:4.5 



is = 7 TeV ! 

LHCb Preliminary - 
— -e c = 0.02 Petern 
— e c = 0.05 Peter.; 
— e c = 0.08 Peter/ 
r c = 0.25 BCFY_ 

a c = 5.0 Kartvel; 

.—.«. = 0.05 CS 




PjL (GeV) 



FIG. 24: Uncertainties in the fragmentation of c 
fragmentation functions from the literature. 
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FIG. 25: Results with overall uncertainties for transverse momentum (left) and rapidity distri- 
butions (right) of Dg for the fcj-factorization approach with the KMR UGDF. For comparison we 
show predictions of other popular approaches. 
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FIG. 26: Transverse momentum distribution of neutral D° mesons for different ranges of rapidities 
specified in the figures. We compare results of the /^-factorization approach with the KMR UGDF 
and those obtained within other approaches known from the literature. 
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FIG. 27: Transverse momentum distribution of charged D + mesons for different ranges of rapidi- 
ties specified in the figures. We compare results of the /ct-factorization approach with the KMR 
UGDF and those obtained within other approaches known from the literature. 
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FIG. 28: Transverse momentum distribution of D* + mesons for different ranges of rapidities 
specified in the figures. We compare results of the /^-factorization approach with the KMR UGDF 
and those obtained within other approaches known from the literature. 
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V. PRODUCTION OF DD PAIRS 



Most of the calculations in the literature concentrates on single meson distributions. We 
wish to focus now on correlation observables for D and D mesons. In order to calculate 
correlation observables for two mesons we follow here, similar as in the single meson case, 
the fragmentation function technique for hadronization process: 

da(pp ->■ DDX) [ D c ^ D (zi) D 5 ^ D (z 2 ) da(pp -> ccX) 



dyidy 2 d 2 p®d 2 p 2t 



I D„ D { Z1 ) D- C ^(z 2 ) da(pp-+ ccX)_ dzidz2 
J zi z 2 dyxdy 2 d 2 p\ t d 2 p\ t 



where: p\ t = p 2 1 = ^ and meson longitudinal fractions Zi,z 2 € (0, 1). The multidimen- 
sional distribution for c quark and c antiquark is convoluted with respective fragmentation 
functions simultaneously. As a result of the hadronization one obtains corresponding two- 
meson multidimensional distribution. In the last step experimental kinematical cuts on 
the distributions can be imposed. Then the resulting distributions can be compared with 
experimental ones. For numerical calculations here we again apply the Peterson model of 
fragmentation function (33| . 

The experimental cross sections for the production of two mesons are also (or even a 
bit more) sensitive to the details of hadronization as it was in the cases of the inclusive 
single D meson production discussed in the previous section. For example in Fig. [29] we 
compare transverse momentum distribution of D° meson provided that D° is also measured 
for two different values of the e c parameter of the Peterson fragmentation function. The 
larger meson transverse momentum, the larger sensitivity to the value of e c . For illustration 
we show the r ang e of transverse momenta relevant for the recent experiments of the LHCb 



collaboration [23] . The effect of the modification of the e c from 0.05 to 0.02 is quite sizeable. 
In the LHCb acceptance, it does not really affect the shape of calculated p t distribution 
but has an important effect for the predictions of integrated cross sections. In Table [TT] we 
compare measured by the LHCb collaboration cross sections for different DD modes with 
our theoretical results. Calculated values for three different UGDFs are consistent with the 
measured ones, taking into account rather large experimental and theoretical uncertainties. 
In particular, this is true only when e c = 0.02 is taken in the calculation of the fragmentation 
process. Results obtained with the KMR UGDF are the closest to the experimental numbers. 



TABLE II: Integrated cross sections for the two mesons modes specified in the table below within 
the LHCb detector. 



Mode fjEXP [nb] KMR t(m c ) Jung set A+ KMS 

e c = 0.05 e c = 0.02 e c = 0.05 e c = 0.02 e c = 0.05 e c = 0.02 



D°D° 


6230 ± 120 ± 630 


r- 1Q q +1346 +654 
oiyo _ g79 _ 576 


6971 


4532 


5814 


2895 


3894 


D°D- 


3990 ± 90 ± 500 


41^ + 1076 + 523 
4100 _ 704 _ 461 


5577 


3626 


4652 


2316 


3115 


D°D S 


1680 ± 110 ± 240 


1/171 + 381 + 185 
14(1 _ 249 _ 163 


1974 


1284 


1647 


820 


1103 


D+D- 


780 ± 40 ± 130 


ooi +215 +105 
ool _ 141 _ 92 


1115 


725 


930 


463 


623 


D+D s 


550 ± 60 ± 90 


coo +152 +74 
000 -99 -65 


790 


513 


659 


328 


441 


D + S D S 




104 + 27 + 13 
1U4 _ 17 _ n 


139 


91 


117 


59 


78 
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FIG. 29: Transverse momentum of D meson within the LHCb acceptance provided that D 
was registered too. Here the KMR UGDF was used. We show results for different values of the 
Peterson fragmentation function parameter e c . 

In Fig. |30]we present transverse momentum distributions of D° meson for the case when 
D°D° pairs are counted. We compare theoretical distributions for different UGDFs (left 
panel) as well as discuss effect of the scale dependence (right panel) on the shape of the p t 
distribution. The experimental data points are normalized by a factor 1/a. The shape of 
the transverse momentum distribution is rather well reproduced by all used UGDFs. The 
normalization, as discussed already in Table [Til is l ess consistent. 




FIG. 30: Uncertainties for the conditional transverse momentum distribution due to the choice of 
UGDF model (left) and due to the choice of scales for the KMR UGDF (right). The experimental 
data of the LHCb collaboration are from Ref . [23| . 

The LHCb collaboration presented also distribution in the D°D° invariant mass M D ojjo. 
In Fig.[3T]we show the corresponding theoretical result for different UGDFs. Both, the KMR 
and KMS UGDFs provide the right shape of the distribution. The dip at small invariant 
masses is due to specific LHCb cuts on kinematical variables. On the other hand the shape 
of the distribution almost does not depend on the choice of the scales for the KMR UGDF 
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(see right panel). 




FIG. 31: Invariant mass distribution of the D°D° system for different UGDFs (left) and uncer- 
tainties due to the choice of the scale for the KMR UGDF (right). The experimental data of the 
LHCb collaboration are taken from Ref. [231 ]. 

The LHCb detector has almost full coverage in azimuthal angle. In Fig. |32] we discuss 
distribution in azimuthal angle between the D° and D° mesons (p D oj=>o. Again the KMR and 
KMS distributions give quite reasonable description of the shape of the measured distribu- 
tion. Both of them, give the enhancement of the cross section at <poo ~ 0. This is due to 
the fact that these approaches include effectively gluon splitting contribution, not included 
in the case of the Jung UGDFs. This was also discussed in Ref. |8j where additional calcula- 
tions of the g*g* — >• gg — >• gcc subprocess in the case of the Jung UGDFs were performed to 
describe azimuthal angle correlation between D and D mesons measured at Tevatron. Some 
dependence of the shape on the choice of the factorization/renormalization scale in the case 
of the KMR UGDF can be also observed (see the right panel). However, still even with the 
KMR UGDF, one can observe some small missing strenght at small angles. It may suggest 
that within the KMR model the gluon splitting contribution is not fully included. 




FIG. 32: Distribution in relative azimuthal angle between D° and D° for different UGDFs (left) 
and uncertainties due to the choice of the scale for the KMR UGDF (right). 
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In order to better understand the result for the azimuthal correlations in Fig. [33] we 
show two-dimensional distributions in invariant mass M D ojjo and azimuthal angle between 
mesons ip D o£>o. The maximum obtained for the KMR UGDF for small relative azimuthal 
angle between D and D mesons corresponds to small invariant masses of the D°D° system. 
This strongly supports the interpretation of the effect as the gluon splitting into cc pair. 
However, one can also see that these interesting shapes of the correlation observable are the 
consequence of the specific LHCb kinematical cuts. 



p p -> D° D° X (s = 7 TeV 




p p -> D° D° X (s = 7 TeV 




p p -> D° D° X (s = 7 TeV 




FIG. 33: Two-dimensional distribution in DD invariant mass and relative azimuthal angle between 
D and D for the KMR and Jung setA+ UGDFs. 



VI. SUMMARY 

First we have discussed the general situation with the cc production at LHC energies. 
We have argued that the cc production is related with small-x physics. Therefore it has a 
good potential to test different models of unintegrated gluon distributions. 

In the present paper we have focused on production of D mesons at the LHC within 
the fc r factorization formalism with unintegrated gluon distributions. Only the Kimber- 
Martin-Ryskin unintegrated gluon distribution gives transverse momentum distributions of 
charmed mesons similar to the recently measured ones by the ATLAS, ALICE and LHCb 
collaborations. Our inclusive theoretical distributions with the KMR UGDFs are very similar 
to those obtained within FONLL or MC@NLO approaches. All other unintegrated gluon 
distributions strongly underpredict the experimental results. This may suggest that some 
mechanism of charm production is still missing. In a followi ng p aper we shall discuss double 



parton scattering effects as a potential missing mechanism [43[ . 

Recently the LHCb collaboration has presented first results for D and D meson two- 
particle distributions. We have presented first theoretical results for such observables. Our 
model calculation with the KMR UGDF relatively well describes both DD meson invariant 
mass distributions as well as DD correlations in relative azimuthal angle between meson and 
antimeson. This shows that the fct-factorization approach is very efficient in describing the 
two-particle distributions. In contrast the NLO QCD approach can be used only in a limited 
region of the phase space but no real results have been presented so far in the literature. 
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